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Abstract 

Dirac fermions on a two-dimensional lattice with disorder are considered. The 
Dirac mass, which controls the gap between the two bands of the fermions, 
is subject to random fluctuations. Another type of disorder is discussed pre- 
sented by a random vector potential. It is shown that the imaginary part 
of the one-particle Green's function can be written as the imaginary part 
of another Green's function which has only poles on the lower half-plane. 
Therefore, it is possible to perform a Cauchy integration for a Lorentzian dis- 
tribution in analogy with the Lloyd model. The results are compared with 
calculations performed in the continuum limit based on renormalization group 

and bosonization methods. 

PACS numbers: 71.20.-b, 71.55.-i, 73.20.Dx 
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I. INTRODUCTION 



The density of states (DOS) in two-dimensional electron systems with a pseudo-gap is 
a subject of interest for a number of physical situations discussed recently . A typical 
model with a pseudo-gap is represented by the Dirac Hamiltonian in two dimensions 

H = zViCTi + iV 2 <72 + m(T S , (1) 

where a a are Pauli matrices including the 2x2 unit matrix ctq- The Dirac equation for 



a state tp then is —dip/dt = Hip. The dispersion relation is E(ki, k 2 ) = ±ym 2 + fcf + k\ 
in the continuum limit. (The lattice will be considered later.) The two signs describe the 
particle and the hole band, respectively. Both bands touch each other if the Dirac mass 
vanishes, as one can see in the DOS p(E) oc \E\Q(E 2 — m 2 ), where O is the step function. 
The touching bands is also a feature of a second order phase transition because the decay 
length of the corresponding Green's function diverges as one goes to the special (critical) 
point m — 0. This behavior is indeed formally related to a number of critical phenomena 
in two-dimensional systems like the ferromagnetic phase transition of the two-dimensional 



Ising model ||10|| . Another physical example, described by the Dirac Hamiltonian, is the 
degenerate semiconductor which exists for m = 0|. Furthermore, the large scale limit of a 
two-dimensional electron gas on a square lattice near the integer quantum Hall transition for 
certain commensurate flux situations (e.g., half a flux quantum per plaquette) is described 
by Dirac fermions [0-0]. A common feature of all these systems is that the DOS at the 
touching bands (i.e., at m — 0) is zero, i.e., there is a pseudo-gap. This raises the question 
whether there is a mechanism which creates states in the pseudo-gap, for instance, electron- 
electron interaction or quenched disorder. This is important in order to understand if there 
is a non-vanishing density of low-energy excitations created by interaction or disorder. In 
this article only the effect of quenched disorder will be analyzed. 

There is a number of studies for the effect of disorder in the pseudo-gap of Dirac fermions. 
A numerical calculation for an electron on a square lattice with half a flux quantum per 
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plaquette shows a non-zero density at low energies || . A coherent potential approximation 
(CPA) of the Dirac fermions with random energy term Ecq added to H also gives a non- 
zero DOS Jl|. A similar result was found for a random mass term ma^ in a modified model 
with N fermion levels per site, using the N oo limit [[11]]. However, these are essentially 



mean-field results which may be affected strongly by fluctuations in the two-dimensional 
system. It is possible that the CPA or N — > oo result are destroyed by fluctuations in 
d = 2. Therefore, as an alternative approach a renormalization group treatment was applied 
to this problem [ID|. From this it turned out that there is asymptotic freedom indicating 



that the pseudo-gap, which is controlled by large scale degrees of freedom, is not affected by 
a random Dirac mass (marginally irrelevant perturbation). However, a rigorous estimation 



leads to a non-zero lower bound of the DOS in the pseudo-gap |12|], at least for a random 
Dirac mass. The renormalization group calculation indicates that the random energy term 
is a relevant perturbation, in agreement with the CPA result. A third type of disorder 
was studied recently by adding a random vector potential to H ||. The renormalization 
group and bosonization treatment indicate a more complicated behavior of the pseudo-gap 
in this case: The average DOS vanishes like \E\ a with a non-universal exponent a > if 
the randomness is weaker than a critical strength. On the other hand, the average DOS 
diverges if the randomness is stronger than the critical strength because of a < 0. (For more 
detailed results see Sect. III.) 

A similar system with a pseudo-gap is the d-wave superconductor. Nersesyan et al. || 
analyzed this system in d = 2 and found for the pseudo-gap of the average DOS p(E) ~ l-^l 1 ^ 7 
in contrast to the linear behavior of the pure system. However, this result is in disagreement 
with others which also find a destruction of the pseudo-gap |J . The effect of disorder in the 
d-wave superconductor will be discussed in a separate article. The aim of this article is to 
present an exact solution for the imaginary part of the single particle Green's function of 
disordered Dirac fermions. 

There are several examples in the theory of a quantum particle in a random potential 
where the average one-particle Green's function can be calculated exactly. Apart from a 



number of one- dimensional examples ||13|| , there is the Lloyd model [p~4| . It is defined by the 
Hamiltonian H = H + V, where H is a Hermitean matrix (e.g., a tight-binding Hamiltonian 
for a particles on a lattice). V is a random potential distributed according to a Lorentzian 
(or Cauchy) distribution 

P(V)dV = (r/n)[r 2 + (V - V Q f]- l dV (2) 

The distribution density has two poles: V = V ±ir. The Green's function G(z) = (H + z)' 1 
must be averaged with respect to the random potential. G(z) as a function of V x for a 
fixed site x is analytic in the upper (lower) complex half-plane if ^sz is positive (negative), 
respectively. Therefore, the path of integration of V x can be closed in that half-plane where 
G is analytic. As a result, only the pole of the Lorentzian density contributes to the integral 
/ GP(V x )dV x . This integration can be performed for all lattice sites leading eventually to 
the average Green's function 

(G{z)) = [H + V + z + isign^z)^- 1 = G(z + isign(Sz)r). (3) 

The average DOS then reads 

(p(E)) = -(l/Nir)\im$?Tr(G(E + ie)), (4) 

e— >0 

where Tr is the trace operator and N is the number of lattice sites. Another examples for an 
exact solution is the DOS of a particle in a homogeneous magnetic field in two dimensions. 
If the corresponding Hilbert space of the particle is projected onto the lowest Landau level 



151 , the average DOS for a white noise potential can be calculated exactly by summing up 
all terms of the perturbation theory with respect to the white noise potential. The exact 
solution is related to the fact that the lowest Landau level system is equivalent to a zero- 
dimensional model. It was discovered by Brezin et al. |16| that the latter is a manifestation 
of the dimensional reduction of the two-dimensional system by 2 due to a supersymmetry of 
the lowest Landau level problem. Unfortunately, the simplicity of the average DOS of the 
lowest Landau level cannot be extended to higher Landau levels. It is also in sharp contrast 



to the complexity of the description of the localization properties [I7|. 



There is some hope that the treatment of an electron on the square lattice in a strong mag- 
netic field is simpler than a continuum model. The lattice model is motivated by numerical 
simulations |18[ and analytic calculations 0-0]. The reason for a simplification is that the 



electron near a quantum Hall transition behaves like a Dirac fermion because the ex- 
citations near the Fermi energy have a linear dispersion. The Hamiltonian of the Dirac 
fermions on a square lattice with unit lattice constant is 

H + iea = (zVi + a)ai + iV 2 o" 2 + mcr 3 + it&o- (5) 

The lattice gradient zV M , with V M /(x) = (1/2) [/(x + e M ) — f(x — e M )] and lattice unit 
vectors e±, e 2 , is Hermitean. Two types of disorder are discussed subsequently: a random 
Dirac mass m and a random vector potential a. The vector potential term is chosen in 
(j|) in the same way as in Ref. It can be considered as a weak disorder approximation 
of a fluctuating Peierls phase factor in Landau gauge. The Green's function now reads 
G(m,ie) = (H + ieuo) -1 , i.e., m and e correspond with the real and imaginary part of z 
in the Green's function of the Lloyd model, respectively. The treatment of this problem is 
rather technical, although the aim is always to find a G with QG = SsG', where the analytic 
properties of G and G' are different: G = (H + ieo~o)~ 1 , as a function of a random variable at 
a given site, has poles on both complex half-planes whereas G' = (H' +ieao)~ 1 has only a pole 
on one of the complex half-plane. The Hamiltonian H' is obtained from H by multiplication 
with a diagonal matrix. The latter depends on the specific type of randomness. 

The article is organized as follows. In Sect. II the random Dirac mass and in Sect. 
Ill a random vector potential are analyzed. The problem of species multiplication due to 
the lattice is discussed in Sect. IV and the projection onto the homogeneous modes on the 
lattice is given. 



II. RANDOM DIRAC MASS 



The matrix H + ieo~o depends on the two complex variables ±m x + it. Thus, in contrast 
to the Lloyd model, G(m, ie) may have singularities in both complex half-planes. Therefore, 
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the Green's function of Dirac fermions is similar to the two-particle Green's function of a 
non-relativistic particle. However, it will be shown subsequently that there is an alternative 
representation for the imaginary part of the Green's function which depends only on a single 
complex variable like the Green's function of the Lloyd model. As a first step, H + ieo~o 
is multiplied by a diagonal matrix Da 3 from the right (D is the staggered diagonal matrix 
D xx > = (— l) xl+X2 5x,a;' with the two-dimensional space coordinates x = (xi,x 2 )) 

H' = i{iVi)Da 2 + i{iV 2 )D(j l + mDa + ieDa 3 , (6) 

where V M -D is Hermitean, since D anticommutes with V^. Hermitean conjugation of H' 
yields 

H ,] = i{iVi)Da 2 + i(iV 2 )Ari + mDa - ieDa 3 . (7) 

Moreover, <t 3 anticommutes with o\ and cr 2 . Consequently, ieDa 3 commutes with all other 
terms in H' . These properties lead to the product 

H'H A = [i{iVi)Da 2 + i(iV 2 )Do- x + mDa ] 2 + eV . (8) 

From the definition of H' follows directly 

H'H' ] = [H + ieor )a 3 Da 3 D{H + iea o y = [H + iea ][H — iea ] = [H — iea ][H + iecr ]. (9) 

The r.h.s. of (§) can also be written H"H"^ with 

H" = i{iV 1 )Da 2 + i{iV 2 )Da t + (mD + ie)a . (10) 

As a result, H" depends only on one complex variable (— l) Xl+X2 m x +ie for a given site x. The 
imaginary part of the Green's function (H +ieao)~ 1 reads (i/2)([H+ieao}~ 1 — [H— ieero] -1 ) = 
e([H — ieao\[H + zecio]) -1 , i.e., it depends on the Hamiltonian only via H 2 . Therefore, the 
identity H 2 + ie 2 a = H"H"^ can be used to write 

~{[H + teo-o]- 1 -[H- tea,]- 1 ) = e([H - iea ][H + ^ao])- 1 = e(H"^H")- 1 

= i -[H"- 1 -{H'«Y\ (11) 
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Thus the imaginary part of the average Green's function (H — ieao)™ 1 can be calculated 
exactly for a Lorentzian distribution due to (|IID , where m x can be integrated out explicitly 
as in the Lloyd model. Only a pole of the distribution contributes leading to the replacements 
e — > t + e and m — > mo in H . This implies 

^([F + Wo]- 1 - [H-iea,]- 1 )) = % -[H^ - {H^\ (12) 

with 

H = zVicri + zV 2 02 + i(e + r)cr + m o"3. (13) 

The imaginary r-term leads always to an exponential decay of the average Green's function 
with a typical decay length £ ~ (m§ + r 2 ) -1 / 2 at e = 0. Moreover, from equ. (||) follows for 
the average DOS 

(p( ie ,m )) = -^Tr[H- 1 ]. (14) 

The dependence on the energy E is obtained from an analytic continuation ie — > ie + E. 
The resulting average DOS is plotted in Fig.l for r = 0.01 and in Fig.2 for r = 0.1. The 
non- vanishing DOS is in agreement with a rigorous proof |12| and a numerical result 0. For 



Gaussian disorder with variance g there is a lower bound Jl2| 

(p(0,0))> Cl e- C2 / 9 . (15) 

with some positive constants c±, c 2 , independent of g. 

In the continuum limit it was argued, using a one-loop renormalization group calculation, 
that random fluctuations of the Dirac mass are irrelevant on large scales |)]||. This implies 
a linearly vanishing DOS and a divergent correlation length at E = m = 0. 

III. RANDOM VECTOR POTENTIAL 

A calculation analogous to that of the random mass system can be performed for a 
random vector potential. For this purpose the orthogonal transformation (o"i + a 3 ) / \pl is 
applied to the Hamiltonian 



H + ieao — > (iVi + a)a 3 — iV2C2 + mai + iea . (16) 

The multiplication of the massless Hamiltonian (i.e., m = 0) from the r.h.s. with D'a^ 
(where D' x x , = {-1) X2 5 X)X ,) yields 

H' = (zVi + a)DV - i(zV 2 )-DVi + ieDV 3 . (17) 



The lattice difference operators iVi-D' and i{i^2)D' are Hermitean. Since -DV3 commutes 
with the first two terms of H', one obtains 

[H + tea ] [H - tea ] = H'H* = [(iVi + a)DV - i(iV 2 )£>Vi] 2 + eV = i/"//" 1 (18) 

with 

H" = (iVi + a)D'<7 - i(iV 2 )I>Vi + +ie(T (19) 

which can be used to establish again equation (|ll|). For (a) = the average imaginary 
part of the Green's function and, therefore, the average DOS for itlq = is related to the 
Hamiltonian H as given in (|13|). 

In contrast to this result, the bosonization of the Dirac fermions in the continuum limit 
leads to a different behavior [H]. For instance, the DOS reads 

(p(E, 0)) ~ £?( 2 -*>/*, (20) 

where z = 1 + Aa/tt (A a is the variance of the fluctuations of the vector potential). I.e., the 
DOS vanishes at E = for z < 2 (weak disorder) and diverges for z > 2 (strong disorder). 
The Green's function behaves like 

(G , x (E,m = 0))~e l ^ x e-^ (21) 

where for E ~ 

A ~ ^-(i-Aa/*) (22) 



6 ~ S-VA^. 
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(23) 



Thus the Green's function decays exponentially for E ^ 0. There is a critical point E = 
where the correlation length diverges with E" 1 . The difference of the results of and the 
present work is probably related to the order of taking the continuum limit and averaging 
over disorder. It is not a consequence of the difference of disorder distributions (Gaussian 
in versus Lorentzian distribution here) because the Gaussian distribution could also be 
treated for H"~-H"i in a strong disorder expansion. The result of this expansion is also 
a finite non-zero DOS and a finite correlation length of G. 



IV. REMARK ON SPECIES MULTIPLICATION 

The phenomenon of species multiplication in a fermion lattice theory is well-known from 
lattice gauge field theories |19| . It is due to several nodes in the energy dispersion of the 



lattice model which indicate the existence of low energy excitations on different length 
scales. The dispersion of the Dirac fermions considered in this article for m = a = is 



E{ki } k,2) = iysin k\ + sin k 2 . It has 9 nodes at kj = 0, ±7r (cf. Fig. 3). In contrast to 



the lattice model the corresponding continuum model, with E(ki,k2) = iy/^i + ^l; has 
low energy excitations only for small wavevectors (i.e., on large scales) as discussed in the 
Introduction. It will be shown in this Section, using the random mass model of Section II, 
that the species multiplication is not the reason for the smooth properties of the one-particle 
Green's function. 

The degeneracy of the low energy behavior of the lattice model can be lifted by introduc- 



ing additional terms in the Hamiltonian JT9|. A possible way is to replace the Hamiltonian H 
by the new Hamiltonian H + S(A — 2)03, where 5 is a positive number (0 < 5 < 1) and A is a 
lattice operator with Af (x) = [f(x+ei)+f(x — ei)+f(x+e2)+f(x — e2)]/2. The dispersion 



of the new H is E(ki, k 2 ) = ±y S 2 (coski + cosfe — 2) 2 + sin k\ + sin k 2 for m = a = 0. 
This is shown in Fig. 4 for 5 = 1/2. It is not clear to the author which tranformation can be 
applied to relate the imaginary part Green's function with the new Hamiltonian in order to 
get the analytic behavior necessary to perform the Cauchy integration with respect to the 
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H 



(24) 



randomness. However, this difficulty can be circumvented by generalizing H to H with 

H + S(A — 2)03 m'o-3 

mV 3 iJ-5(A + 2)a 3 , 

where m' is a random variable which is statistically independent of m with mean zero. Now 
the orthogonal transformation 



1 / Co 00 
^ \0o -0o. 



(25) 



(26) 



rotates the diagonal part (Acr 3 , — Acr 3 ) in the off-diagonal positions and the off-diagonal part 
into the diagonal position (mV 3 , — mV 3 ) such that 

H - (25 - m> 3 Acx 3 

Ao- 3 H - (25 + m')a ?J , 

The random variables M x = —25 + m x + m' x and M^, = —25 + m a . — m^, in the diagonal part 
of H can now be considered as new independent random variables. 
The transformation 

/o-o \ (Da 3 

\0 -00/ V Da 3 
HDa 3 - (25 - m')Da ADa 



-HDa 3 - (25 + m')Da 0i 



H' 



DAa 

generates the Hermitean matrix H' . 

Using for the r.h.s. of (|27D the notation T Q HTi and applying the property 

ToHT, = TiHT 



(27) 



(28) 



one obtains 

(H' - ieD l3 )(H' + ieD 73 ) = T (H - ze 7 o)T 1 T (F + ie lo )T x = T (H - ie lo )(H + ze 7o )T 

(29) 

with the diagonal matrices 70 = (cro,o"o) and 7 3 = (cx 3 , — cr 3 ). Moreover, one has for the 
imaginary part of the one-particle Green's function as before (i/2) (H + ie7o) _1 — (H — 



^7oJ 



(H — ie^o)(H + ie7o) . Due to (E3) and T 1 = T this can be rewritten as 
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eT \(H'-teD l3 )(H' + teD l3 )} T . 



(30) 



The l.h.s. of (M) reads 



[H' - ieD l3 )(H' + ieD l3 ) = {H'f + e 2 7 o = (#' ~ iej ){H' + *e 7o ) 



(31) 



because H' and D^ 3 commute. This implies for ( |30| 



T = - 2l 



Consequently, the imaginary part of the Green's function satisfies 



(H + ie-fo)- 1 - (H-ie-f 



(^ / + ze 7o )- 1 -(^ / -ze 7 or 1 ]T , 



(32) 



(33) 



analogously to (|TTD . At a given site x the matrix H'+ie'yo depends on the random variables in 
the combinations (— l) Xl+X2 M x +ie and (— l) Xl+X2 M' x +ie. Assuming a Lorentzian distribution 
for M x and M' x , the integration can be performed again as in Sect. II. As a result the 
imaginary part of the average Green's function is 

H + 5(A-2)a 3 N ~' 

H - 5(A + 2)a 3 



(34) 



where H is the average Hamiltonian (0). Thus the lifting of the degeneracy of the nodes 
in the dispersion relation does not change the analytic behavior of the average one-particle 
Green's function. 



V. CONCLUSION 

An exact expression for the average imaginary part of the one-particle Green's function 
and the average DOS of two-dimensional lattice Dirac fermions have been derived for a ran- 
dom Dirac mass and for a random vector potential. We have shown that there is a non-zero 
DOS due to disorder and there is a finite decay length for the average one-particle Green's 
function. This implies the creation of a non- vanishing density of low-energy excitations due 
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to disorder in a vicinity of E = M = 0. These lattice results are in agreement with numeri- 
cal simulation [[3]. However, they are in disagreement with the results of a renormalization 
group calculation and a bosonization approach for a continuous system of Dirac fermions 
PH, where the DOS vanishes or diverges at E — M — 0. Moreover, the lattice model 
does not exhibit the critical properties of the Green's function and the DOS found in the 
renormalization group calculation and in the bosonization approach. It is possible to take 
the continuum limit of the lattice model after the averaging over disorder, for instance, in 
the Hamiltonian flT3|). This, however, does not lead to a critical behavior. It seems that 
the critical behavior of the DOS is a consequence of taking the continuum limit first and 
performing the averaging over disorder afterwards. This is plausible because the effect of 
randomness is much stronger in the continuum due to statistically independent fluctuations 
on arbitrarily short scales. It is shown in Sect. IV that species multiplication, which is a 
special effect of the lattice model, is not the reason for the smooth behavior of the average 
DOS. 
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FIGURES 



FIG. 1. Average density of states for disorder strength r = 0.01. 

FIG. 2. Average density of states for disorder strength r = 0.1. 

FIG. 3. Energy dispersion E(k±,k2) of lattice Dirac fermions. 
FIG. 4. Energy dispersion of lattice Dirac fermions after lifting the degeneracy of the low energy 

properties 
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Fig.2 
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Fig.3 




Fig.4 



